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Abstract — Quantum convolutional codes can be used to protect 
a sequence of qubits of arbitrary length against decoherence. 
We introduce two new families of quantum convolutional codes. 
Our construction is based on an algebraic method which allows 
to construct classical convolutional codes from block codes, 
in particular BCH codes. These codes have the property that 
they contain their Euclidean, respectively Hermitian, dual codes. 
Hence, they can be used to define quantum convolutional codes 
by the stabilizer code construction. We compute BCH-like bounds 
on the free distances which can be controlled as in the case of 
block codes, and establish that the codes have non-catastrophic 
encoders. 

I. Introduction 

Quantum convolutional codes provide an alternative to 
quantum block codes to protect quantum information for 
reliable quantum communication. Ollivier and Tillich launched 
the stabilizer framework for quantum convolutional codes [11]. 
Using this stabilizer framework Forney et al. constructed rate 
(n— 2) / n quantum convolutional codes [3]. Recently, two of us 
constructed quantum convolutional codes from product codes 
[5] and derived an algorithm to construct non-catastrophic 
encoders and encoder inverses [6]. In [1], a generalized 
Singleton bound for a class of quantum convolutional codes 
has been established, together with a family of codes based 
on generalized Reed-Solomon codes meeting this bound. 

Unit memory convolutional codes are an important class 
of codes that appeared in a paper by Lee [10]. He also 
showed that these codes have large free distance df among 
other codes (multi-memory) with the same rate. Convolutional 
codes are often designed heuristically. However, classes of 
unit memory codes were constructed algebraically by Piret 
based on Reed-Solomon codes [12] and by Hole based on 
BCH codes [8]. In a recent paper, doubly-cyclic convolutional 
codes are investigated which include codes derived from Reed- 
Solomon and BCH codes [4]. These codes are related, but not 
identical to the codes defined in this paper. 

The main results of this paper are: (a) a method to construct 
convolutional codes from block codes (b) a new class of 
convolutional stabilizer codes based on BCH codes. These 
codes have non-catastrophic dual encoders making it possible 
to derive non-catastrophic encoders for the quantum convolu- 
tional codes. 

II. Background 

A. Convolutional Codes 

We briefly recall the basic facts about classical convolutional 
codes relevant for our discussion. Let F g be a finite field with 



q elements. A convolutional code C of length n and dimension 
k over F g is a free module of rank k that is a direct summand 
of F q [D} n . A matrix G in F q [D] kxn such that C = imG = 
{uG | u G F g [Z3] fe } is called a basic generator matrix of G, 
and a matrix H G F 9 [£>]("- fc ) xn such that C = keriJ* = 
{v | v G F q [D] n , v_ff* = 0} is called a basic parity check 
matrix of C. 

The existence of a convolutional code C is equivalent 
to the existence of four matrices G G F q [D] kxn , H G 
F q [D}( n ~V xn , K G F q [D] nxk , and L G F q [D] nx( * n ~V such 
that C = imG = keriJ*, GK = l Fg[D p, and iSH 1 = 

lF,[D]"- fc = HL. 

Let Vi denote the maximum of the degrees among the 
polynomials in the ith row of a basic generator matrix G, 
and let the memory m be the maximal value of v.i. A basic 
generator matrix of a convolutional code C is called reduced 
if the overall constraint length v = V\ + • • • + Uh has the 
smallest value among all basic generator matrices of C. It 
is often convenient to express the generator matrix as G = 
G + GiL> + • • ■ + G m D m , where G 4 G F k q xn . 

Let F q ((D)) be the field of Laurent series consisting of 
elements of the form v(D) — J2i y i-D l for Vi G F q and 
Vi = for i < r for some r G Z. We associate with a con- 
volutional code G another module G°° = {u(D)G | u(D) G 
F q ((D)) k }, The entries of a generator matrix G of C°° might 
be rational functions. Let v(D) = (vx(D), . . . , v n (D)) G 
F q ((D)) n where Vi(D) = ~}2,jVijD^. Then we can identify 
v(D) with an element in F™((D)) as ^ ■ VjD J , where Vj = 
(vij , . . . , v n j) G F™. We define the weight of v(D) as 
wt(v(Z?)) = X)iez w K Vi )- A generator matrix G is called 
catastrophic if there exists a u(D) G F q ((D)) k of infinite 
Hamming weight such that u(D)G G G°° has finite Hamming 
weight. The free distance df of C is defined as 

d f = rnin{wt(v(L>)) | v(D) G C,v(D) jt 0}. (1) 

A rate k/n convolutional code with memory m, overall 
constraint length v, and free distance df is denoted by 
(n,k,iy;m,df) q . Sometimes a shorter notation (n 1 k,v) q is 
also used. 

The Euclidean inner product of two n-tuples u(D) = 
^i^D' and V (-C) = T,j v tD j in F q [D] n is defined as 
(u(D)\v(D)) — u i ' v i- The Euclidean dual of a con- 
volutional code C is denoted by G x = {u(D) G F q [D] n \ 
(u(D)\v(D)) = for all v(D) e C}. Note that H(D), the 
parity check matrix of C, does not generate the Euclidean 
dual in general. Instead, one has to reverse the order of 



the coefficients of the polynomials in H(D), i.e. consider 
the matrix D m H(l/D), where m is the memory of the 
code generated by H(D). For codes over F g 2, we define the 
Hermitian inner product as (u(D)\v(D))h — J2i u i' v h wnere 
vu,Vi GF™ 2 and vf = «,.. 
C is then C ±h = {u(D) G 
for all v(D) G C}. 



J. The Hermitian dual of 

[Df | (up)|v(U)) h - 



B. Quantum Convolutional Codes 

We briefly describe the stabilizer framework for quantum 
convolutional codes, see also [1], [7], [11]. The stabilizer is 
given by a matrix 



S(D) = (X(D)\Z(D))eF q [D] 



(n — k)x2n 



(2) 



which satisfies the symplectic orthogonality condition = 
X(D)Z(1/ 'Df - Z(D)X(1/ '£>)*. Let C be a quantum convo- 
lutional code defined by a stabilizer matrix as in eq. 10. Then 
n is called the frame size, k the number of logical qudits per 
frame, and k/n the rate of C. It can be used to encode a 
sequence of blocks with k qudits in each block (that is, each 
element in the sequence consists of k quantum systems each 
of which is q-dimensional) into a sequence of blocks with n 
qudits. 

The memory of the quantum convolutional code is defined 
as m — maxi<i< n -k,i<j< n (max(degXij(D),degZij(D))). 
We use the notation \(n,k,m)] q to denote a quantum con- 
volutional code with the above parameters. We can identify 
S(D) with the generator matrix of a self-orthogonal classical 
convolutional code over F q or F q 2, which gives us a means 
to construct convolutional stabilizer codes. Analogous to the 
classical codes we can define the free distance, df and the 
degree v, prompting an extended notation \{n, k, m; v, df)] q . 
All the parameters of the quantum convolutional code can 
be related to the associated classical code as the following 
propositions will show. For proof and further details see [ 1 L 1 ! 

Proposition 1: Let (n, (n — k)/2, v; m) q be a convolutional 
code such that C < C , where the dimension of C is 
given by (n + k)/2. Then an [(n, k, m; v, df)\ q convolutional 
stabilizer code exists whose free distance is given by df — 
wt(C \C), which is said to be pure if df = wt(C ). 

Proposition 2: Let C be an (n,(n — k)/2,i>;m) q 2 con- 
volutional code such that C C C ±h . Then there exists an 
[(n, k, m; u, dt)] q convolutional stabilizer code, where df = 
wt(C- L " \ C). 

III. A Construction of Convolutional Codes 

In this section, we give a method to construct convolutional 
codes from block codes. This generalizes an earlier construc- 
tion by Piret [13] to construct convolutional codes from block 
codes. One benefit of this method is that we can easily bound 
the free distance using the techniques for block codes. Another 
benefit is that we can derive non-catastrophic encoders. 

'A small difference exists between the notion of memory defined here and 
the one used in [1]. 



A. Convolutional Codes from Block Codes 

Given an [n,k,d] q block code with parity check matrix H, it 
is possible to split the matrix H into m+1 disjoint submatrices 
Hi, each of which has n columns, such that 



H = 



Ho 
Hi 

H m 



Then we can form the polynomial matrix 
G(D) =H a + HiD + H 2 D 2 + . . . + 



H m D r 



(3) 



(4) 



where the number of rows of G{D) equals the maximal 
number k of rows among the matrices H. The matrices Hi 
are obtained from the matrices H by adding zero-rows at the 
bottom such that the matrix Hi has k rows in total. Then 
G(D) generates a convolutional code. The fact that the Hi 
come from a common block code allows us to characterize 
the parameters of the convolutional code and its dual using 
the techniques of block codes. Our first result concerns a non- 
catastrophic encoder for the code generated by G(D). 

Theorem 3: Let C C be an [n, k, d] q linear code with 
parity check matrix H G Fq™ fc ) x ™. Assume that H is 
partitioned into submatrices Ho, Hi, ... , H m as in eq. (0 such 
that k = rkffo and rk Hi < k for 1 < i < m. Define the 
polynomial matrix G(D) as in eq. ©. Then we have: 

(a) The matrix G(D) is a reduced basic generator matrix. 

(b) If the code C contains its Euclidean dual C , respectively 
its Hermitian dual C ±h , then the convolutional code V — 
{v(L>) = u(D)G(D) | u(D) G F^- k [D}} is contained 
in its dual V , respectively its Hermitian dual V ±h . 

(c) Let df and dj- respectively denote the free distances of V 
and V . Let di be the minimum distance of the code Ci = 
{v £ FJ vHf = 0}, and let d 1 - denote the minimum 
distance of C . Then the free distances are bounded by 
min{<io + d m , d} < dy 1 < d and df > d ± . 

Proof: To prove the claim (a), it suffices to show that (i) 
G(0) has full rank k, (ii) (coeff(G(D) i3 , £> I/ *))i<i<«,i<j<n- 
has full rank k, where for f(D) = 2~2i>o a iD l we define 
cocff (f(D), D l ) — di, and (iii) G(D) is non-catastrophic; 
cf. [12, Theorem 2.16 and Theorem 2.24]. 

By definition, G(0) = Ho has rank k, so (i) is satisfied. 
Condition (ii) is satisfied, since the rows of H are linearly 
independent; thus, the rows of the highest degree coefficient 
matrix are independent as well. 

It remains to prove (iii). Seeking a contradiction, we assume 
that the generator matrix G(D) is catastrophic. Then there 
exists an input sequence u(D) — X^ 11 *-^ 1 e Fg((-C)) K 
with infinite Hamming weight that is mapped to an output 
sequence v(D) = u(D)G = J2 t v l D z £ F q ({D)) n with finite 
Hamming weight, i.e. Vj = for all i > iq. We have 



U 



i+mHo 



U 



i+m 



_ X H 



Uji/ m , (5) 



where v J+m 6 F" and Uj G F£. By construction, the vector 
spaces generated by the rows of the matrices Hi intersect 



trivially. Hence Vj = for i > iq implies that Ui-jHj = 
for j = 0, . . . , to. The matrix Hq has full rank. This implies 
that Ui = for i > io, contradicting the fact that u(D) has 
infinite Hamming weight; thus, the claim (a) holds. 

To prove the claim (b), let v(D) = J^i v iD l , w(D) = 
J^WiD 1 be any two codewords in V C F"[Z?]. Then from 
eq. (O, we see that v* and Wj are in the rowspan of H i.e. 
they are elements of C , for any i,j G Z. Since C CC = 
(C ) , it follows that Vj • w 3 = 0, for any i,j G Z which 
implies that (v(£>)|w(£>)) = X),: e z v » ' w « = °- Hence V C 
7 X . Similarly, we can show that if C^ h C C, then C V^ h . 

For the claim (c), without loss of generality assume that the 
codeword c(D) = Yli=o c iD l is in V 1 -, with Co ^ ^ q . 

It follows that {D' l c{D)\D l G 3 {D)) = for i, Z > 0, where 
Gj(D) denotes the jth row of G(D). In particular we have 
co-ff^j = and c^7?q = 0. It follows that Co G C m and eg G 
Co. If ^ > 0, then wt(co) > d m and wt(c^) > do implying 



wt(c(L>)) > d + d m . If £ = 0, then (Z>*co|Gj(D)) 







implies co-ff' = for < i < TO, whence coi?* = and 
Co G C, implying that wt(co) > d. It follows that wt(c(D)) > 
min{c?o + d m , d}, giving the lower bound on dj-. 

For the upper bound note that if Co is a codeword of C, 
then coH\ = 0. From Co we can construct a codeword c(D) 
by padding with zeros. Now, {D l c(D)\D l G 3 {D)) = and 
hence c{D) G V^. Since wt(c(D)) = wt(co) we obtain that 
dj < d. 

Finally, let c(D) = £\ c iD l be a non-zero codeword in V. 
We saw earlier in the proof of (b) that every c» is in C . Thus 
d f > min{wt(ci) | a ^ 0} > d ± . ■ 

A special case of our claim (a) has been established by a 
different method in [8, Proposition 1]. 

IV. CONVOLUTIONAL BCH CODES 

One of the attractive features of BCH codes is that they 
allow us to design codes with desired distance. There have 
been prior approaches to construct convolutional BCH codes, 
see [8], [14], and most notably [4], where one can control 
the free distance of the convolutional code. Here we focus 
on codes with unit memory. Our codes have better distance 
parameters as compared to Hole's construction [8] and are 
easier to construct compared to [14]. 

A. Unit Memory Convolutional BCH Codes 

Let F g be a finite field with q elements, n be a positive 
integer such that gcd(n, q) = 1. Let a be a primitive nth root 
of unity. A BCH code C of designed distance 5 and length n is 
a cyclic code with generator polynomial g(x) in F q [x]/(x n — 
1) whose defining set is given by Z = CbUCb + iU- • -UCb+s-2, 
where C x — {xq l mod n | i G Z, i > 0}. Let 

„6 ^,26 „,b(n-l) 



S,b 



Then C 



a 

•V&+1 



a 

V 2(M 



a b{n 
Jb+X){n-1) 



1 a (b+& -V a 2(b+5 - 2) 



a 



(6+5-2)(n-l) 



{v G F' Q 



0}. If 



ord n (g), 



every entry in the matrix Hs b as a column vector over some 
F 9 -basis of F q r, and removing any dependent rows. Let 
B = {bi, . . . ,b r } denote a basis of F g r over F g . Suppose 
that w = (wx, . . . ,w n ) is a vector in F^r, then we can 
write Wj = + • • • + Wj :T b r for 1 < j < n. Let 

w"^ = . . • , w n .i) be vectors in F™ with 1 < i < r, 

For a vector v in F£, we have v • w = if and only if 
v • W M = for all 1 < i < r. 

For a matrix M over F 9 r, let exs(M) denote the matrix 
that is obtained by expanding each row into r rows over F q 
with respect to the basis B, and deleting all but the first rows 
that generate the rowspan of the expanded matrix. Then H = 
ex B (H St b)- 

It is well known that the minimum distance of a BCH code 
is greater than or equal to its designed distance S, which is 
very useful in constructing codes [9]. Before we can construct 
convolutional BCH codes we need the following result on the 
distance of cyclic codes. 

Lemma 4: Let gcd(n, q) — 1 and 2 < a < (3 < n. Let 
C C F^ 1 be a cyclic code with defining set 

Z = {z | z G C Xl a < x < /3, x ^ mod q}. (6) 
The minimum distance A(a,(3) of C is lower bounded as 



A(a,/3) > 



q + |X8 - a + 3)/q\ - 2, if p - a > 2q - 3; 
L(/3 - a + 3)/2j , otherwise. 



Proof: Our goal is to bound the distance of C using the 
Hartmann-Tzeng bound (for instance, see [9]). Suppose that 
there exists a such that A = {z, z + 1, . . . , z + a — 2} C Z. 
Suppose further, that there exists b, where gcd(6, q) < a and 
A + jb = {z + jb, z + 1 + jb, . . . , z + a - 2 + jb} C Z for 
all < j < s. Then by [9, Theorem 4.5.6], the minimum 
distance of C is A (a, (3) > a + s. 

We choose b = q, so that gcd(n, q) = 1 < a is satisfied for 
any a > 1. Next we choose A C Z such that \A\ = q — 1 
and A + jb C Z for < j < s, with s as large as possible. 
Now two cases can arise. If (3 — a + 1 < 2q — 2, then there 
may not always exist a set A such that \A\ = q — 1. In this 
case we relax the constraint that \A\ = q — 1 and choose A 
as the set of maximum number of consecutive elements. Then 
\A\=a-l> [{(3 - a + 1)/2J and s > giving the distance 
A(a, (3)>[{(3-a + 1)/2J + l=\fi-a + 3)/2j. 

If [f3 — a + 1) > 2q — 2, then we can always choose a 
set A G {z | a < z < a + 2q — 3, z ^ mod (?} such that 
|yl| = q — 1. As we want to make s as large as possible, the 
worst case arises when A = {a + q — 1, . . . , a+2q — 3}. Since 
C Z holds for < j < s, it follows a+2g-3 + sq < j3. 
Thus s < [(/? - a + 3)/q\ - 2. Thus the distance A(a, /3) > 
9 + L(a - + 3)/q\ - 2. M 

Theorem 5 (Convolutional BCH codes): Let n be a posi- 
tive integer such that gcd(n, q) = 1, r = ord n (q) and 
2 < 25 < <5 max , where 



then a parity check matrix H for C is given by writing 



"(9 



k/21 



l-( ? -2)[r odd]) 



Then there exists a unit memory rate k/n convolutional BCH 
code with free distance df > 5+l+A(5+l, 25) and k = n—n, 
where n — r [5(1 — 1 / . The free distance of the dual is 

> tax + 1. 

Proof: Let C C F™ be a narrow-sense BCH code of 
designed distance 25 + 1 and let B a basis of F qr over F q . 
Recall that a parity check matrix for C is given by H = 
exs(i?25+i,i)- Further, let H = exs(Hs+i } i), then from 



H 



2<5+l,l 



-ff ,5+1, 5+1 



(7) 



it follows that H = [Hq,H{] , where Hi is obtained 
from exs{H$ + i ,5+1) by removing all rows common to 
exB(Hs+i.i). The code Co with parity check matrix Hq = 
exs(iJ ( 5 + i ! i) coincides with the narrow-sense BCH code of 
length n and designed distance 5 + 1. 

By [2, Theorem 10], we have dimC = n — r [2(5(1 — l/q)~\ 
and dim Co = n — r [5(1 — 1/q)] which implies rkH = 
r [25(1 - l/g)~|, rkiJo = r [5(1 - 1/q)], and rkffi = 
rkH~rkH a = r [2(5(1 - 1/q)] -r [5(1 - 1/q)]. For a; > 0, 
we have \x] > \2x] — \x]; therefore, k — rk-ffo > rki?i. 

By Theorem OaJ, the matrix _ff defines a reduced basic 
generator matrix 



G(D) = H Q + DH X 



(8) 



of a convolutional code of dimension k, while its dual which 
we refer to as a convolutional BCH code is of dimension n— re. 

Now Hi is the parity check matrix of a cyclic code, C\ 
of the form given in Lemma [4] i.e. the defining set of C\ is 
Z\ as defined in © with a = 5 + 1 and /3 = 25. Since Hi 
is linearly independent of Hq we have x ^ mod (7 in the 
definition of Z\. 

By Theorem [3jlcJ, the free distance of the convolutional 
BCH code is bounded as min{do + di,d} < df < d. 
By Lemma g] di > A(5 + 1, 25) and by the BCH bound 
do > 5 + 1. Thus d f > 8 + 1 + A(5 + 1, 25). The dual free 
distance also follows from Theorem [5101 as dj- > d . But 
d 1 - > 5 max + 1 by [2, Lemma 12]. ■ 



V. Constructing Quantum Convolutional Codes 

Under some restrictions on the designed free distance, we 
can use convolutional codes derived in the previous section to 
construct quantum convolutional codes. 

Theorem 6: Assume the same notation as in Theorem [5] 
Then there exists a quantum convolutional code C with pa- 
rameters [(n, n — 2k, l)] q , where n = r [5(1 — 1/q)]. For the 
free distance of C the bound df > 5 + 1 + A(5 + 1, 25) holds 
and it is pure to d! > 5 max + 1. 

Proof: We construct a unit memory (n, n — n) q classical 
convolutional BCH code as per Theorem [5] Its polynomial 
parity check matrix G(D) is as given in eq. ([§}. Using the 
notation as in the proof of Theorem [5] we see that the 
code contains its dual if H is self-orthogonal. But given 
the restrictions on the designed distance, we know from [2, 
Theorem 3] that the BCH block code defined by H contains its 
dual. It follows from Theorem[3j|b]l that the convolutional BCH 



code contains its dual. From Proposition 03 we can conclude 
that there exists a convolutional code with the parameters 
[{n,n — 2k, l)] g . By Theorem [5] the free distance of the dual 
is d! > 5 max + 1, also implying its purity. ■ 

Another useful method to construct quantum codes makes 
use of codes over F q 2 . 

Theorem 7: Let 2 < 25 < [n{q r - l)/(q 2r - 1)J, where 
and r = ord„(q 2 ). Then there exist quantum convolutional 
codes with parameters [(n, n — 2k, l)] q and free distance df > 
5 + 1 + A(5 + 1,25), where re = r [5(1 - 1/q 2 )] . 

Proof: By Theorem [5] there exists an (n, n — k, l) g 2 
convolutional BCH code with the polynomial parity check 
matrix as in eq. dHJ. The parent BCH code has design distance 
25+1 and given the range of 5, we know by [2, Theorem 14] 
that it contains its Hermitian dual. By Theorem [Mb}, the 
convolutional code also contains its Hermitian dual. By Propo- 
sition [2] we can conclude that there exists an [(n, n — 2k, l)] q 
code with df > 5 + 1 + A(5 + 1,25). ■ 
We conclude by noting that the convolutional codes in The- 
orems [6] and [7] have non-catastrophic encoders and encoder 
inverses. This follows directly from the fact that G(D) in 
eq. (O is a basic generator matrix (cf. [6], [7]). 
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